Stark tuning of the charge states of a two-donor molecule in silicon 
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Gate control of phosphorus donor based charge qubits in Si is investigated using a tight-binding 
approach. Excited molecular states of P2+ are found to impose limits on the allowed donor sepa- 
rations and operating gate voltages. The effects of surface (S) and barrier (B) gates are analyzed 
in various voltage regimes with respect to the quantum confined states of the whole device. Effects 
such as interface ionization, saturation of the tunnel coupling, sensitivity to donor and gate place- 
ment are also studied. It is found that realistic gate control is smooth for any donor separation, 
although at certain donor orientations the S and B gates may get switched in functionality. This 
paper outlines and analyzes the various issues that are of importance in practical control of such 
donor molecular systems. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 



I. INTRODUCTION 

Single donor systems in silicon have been the subject 
of much research in recent years from both theorists and 
experimentalists alike. While a number of experiments 
have probed into the physics of single donors [H,I3jBSj@1j 
many others have concentrated on precision control and 
fabrication of qubits based on individual donors @, 0|- 
Apart from the initially predominant effective mass based 
approaches for donor modeling [H, 0] , a number of other 
techniques have been successfully applied to study real- 
istic devices and interpret experimental results in detail 
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Although the initial qubit proposals utilized the s pin 
of the donor electron or nucleus to encode qubits HI, 14 , 
1 5( | , one can also envisage fabricating and controlling the 
charge degrees of freedom of a simple donor molecule: a 
singly ionized double donor system (ljj . Steps along this 
direction have been achieved with charge state transfer 
between two P donors in silicon reported [1 71 ] , and a sin- 
gle electron in a donor-interface double well forming a 
new hybrid molecular system Despite possessing the 
disadvantages of shorter decoherence times, such a charge 
qubit system is easier to probe experimentally and may 
help to demonstrate the feasibility of theoretical quantum 
computation concepts. The donor based charge qubit is 
also more amenable to measurements as the donor elec- 
tron can be localized to particular impurities with rel- 
ative ease by means of appropriate gate placement and 
voltage pulses. Microwave driven experiments have al- 
ready been proposed to investigate the parameter space 
and physical operation of such devices [l8|. Precision 
placement of a few donors has also been achieved re- 
cently, highlighting the urgency for a detailed theoretical 
study of a full charge qubit device before its physical re- 
alization is achieved [17J] . Charge qubits also form the 
building blocks of the coherent charge transfer mecha- 



nism proposed in Ref [191 ]. In such a scheme, quantum in 
formation encoded in the spin or in the charge of a donor 
electron can be transported coherently through a chain of 
ionized donors by an adiabatic pathway realized through 
voltage pulses applied to electrical gates. Termed as Co- 
herent Tunneling Adiabatic Passage (CTAP), this mech- 
anism may provide a robust way to transfer information 
in a circuit of donor spin qubits [13, HH ■ 

An important parameter of interest in a P2+ molecule 
is the tunnel coupling between the two impurity quan- 
tum dots defined by the coulomb potential of each donor 
nucleus. This coupling is expressed as the difference be- 
tween the two lowest eigen states of the molecular sys- 
tem. A previous work [221 ] had utilized Kohn-Luttinger 
p3| type effective mass based variational envelope func- 
tions modulated by Bloch states of the six conduction 
band valleys of Si to show that the tunnel coupling suffers 
from the same sensitivity to relative donor placements as 
the inter-donor exchange coupling for spin qubits. In Ref 
j2~il ] , the same approach was used to show that the tunnel 
coupling can be controlled smoothly by constant electric 
fields for both homo and hetero polar donor species. As 
experiments are getting closer to fabricating an actual 
donor based charge qubit system 17 1. it is important to 
develop a detailed theoretical model that analyzes the 
various aspects of qubit control and design space, and 
predicts accurate numbers for quantities likely to be mea- 
sured by experiments. Our goal in this paper is to serve 
this very purpose by modeling the P2+ molecule with a 
comprehensive numerical approach. 

We employ atomistic tight-binding theory with a 
model for P impurity in Si that accurately captures the 
basic single donor physics such as the valley-orbit split 
donor spectrum. By solving the full tight-binding Hamil- 
tonian for realistic systems of about 3 million atoms, we 
are able to obtain any number of states of the P2+ spec- 
trum. Previous works 22j, [24j on donor charge qubits 
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have ignored the effects of the excited manifold on device 
operations. However, in the presence of significant gate 
bias, these excited states can disrupt normal device oper- 
ation by entering the manifold of states used in quantum 
computing operations. Furthermore, effects of realistic 
gate potential profiles on the tunnel coupling need to be 
investigated as opposed to the simplistic constant elec- 
tric fields used in earlier works. A gate can also cause 
surface ionization of the bound electron at higher biases, 
thus limiting the control regime in practice. Presence 
of nearby interfaces can significantly distort the donor 
wave functions and affect charge qubit operations. Since 
the previous works concentrated on bulk systems, these 
effects have been neglected. 

In our tight-binding model coupled with electrostatic 
gate potentials from a commercial Poisson solver [251 ] , we 
can investigate all these effects under one framework with 
considerable accuracy. It is also important to go beyond 
the effective mass theory (EMT) assumptions of only 
the valley minima states contributing to the donor wave 
functions, and to consider a more comprehensive Bloch 
structure of Si as is done in tight-binding. In the case 
of inter-qubit exchange coupling for spin qubits, Wellard 
[2^ showed that oscillations of the exchange coupling 
J(V) could be damped to some extent if an extended 
set of Bloch states is considered rather than the six val- 
ley minima states to expand the donor wave function. 
Wellard ct al. [2(| also calculated the angular depen- 
dence of J(V) for a fixed radial separation of donors, and 
showed that a gate bias was not able to alter J(V) signif- 
icantly for certain angular separations between donors. 
It remains to be seen if the tunnel coupling also suffers 
from this controllability problem. 

This paper is organized into three parts. In Section 
I, the geometry of a P donor based charge qubit in Si 
is described, and the various control parameters studied 
here are outlined. In Section II, the details of the tight- 
binding method is described. Section III presents the 
results and discusses the controllability issues in detail. 
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FIG. 1: The schematic of a donor based charge qubit device 
showing the various design parameters. The surface gate volt- 
age Vs detunes the charge states of the donor molecule, while 
the barrier gate voltage Vb controls the tunneling barrier be- 
tween the donors. 



Dl from the lateral interfaces. Typical devices have di- 
mensions 50 nm x 40 nm x 30 nm. While a larger device 
volume would be ideal, the compute times increase much 
more making the problem intractable. This device vol- 
ume, however, captures most of the Physics of the double 
donor molecule. 

First, we study the molecular spectrum of the system 
without any gate voltage taking into account the compli- 
cated band structure of Si. Then we investigate how this 
spectrum, and in particular the tunnel coupling, can be 
controlled by surface and barrier gates. Effects of design 
parameters such as Ls, Dy, R, 6, Vs, and Vb are also 
explored. In particular, the energy gaps A12 and A23 are 
studied. The energy gap A 12 represents the energy dif- 
ference between the 1st excited state (S2) and the ground 
state (Si), while A23 represents the energy difference be- 
tween the 2nd excited state(Ss) and the 1st. The device 
schematic and notations of Fig 1 is used throughout this 
paper. The molecular states arc labeled as Si, S2, S3, 
and so on, from more strongly bound to less (Si being 
the ground state). 



II. P 2 + DONOR DEVICE GEOMETRY 



III. TIGHT BINDING SOLUTION OF THE 
SOLID-STATE DONOR MOLECULAR SYSTEM 



A schematic of a P2+ donor molecular device in [001] 
grown Si is illustrated in Fig 1. The P donors are sepa- 
rated by a distance R nm and an angle O measured from 
the [100] axis. The donors form a molecular system anal- 
ogous to H 2 + except the Si band structure complicates 
the scenario. A barrier gate (B-gate) is placed midway 
between the impurities and controls the potential bar- 
rier between them. A surface gate (S-gate) is placed a 
distance L5 (measured with respect to the center of the 
S-gate) away from the left impurity, and controls the de- 
tuning of the impurity states, the oxide thickness T ox is 
5 nm, while the gate lengths arc 10 nm. The impurities 
arc buried at a depth Dy below the oxide and a distance 



The approach used here is the semi-empirical tight- 
binding (TB) method [13] with the 20 band sp 3 d 5 s> 

In this method, the Hamilto 



nearest-neighbor model, 
nian is expressed in real space with a basis of localized 
atomic orbitals. The 20 band TB parameters for Si are 
optimized by genetic algorithm [28|, [2^] to reproduce the 
bulk band structure of the host. Once a set of such pa- 
rameters are found, it can be used for atomistic modeling 
of a generic device made of the host. A single P impurity 
is represented by a coulomb potential screened by the 
dielectric constant of Si (est) an d subjected to a cut-off 
potential Uq at the impurity site. This core correcting 
potential Uq is adjusted to reproduce the single donor 
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ground state binding energy. The total TB Hamiltonian 
is of the form. 



H = H Q + eVbi + eV D2 + - — — 

47T€Si|-Rl - i?2 

+ eV G (V s ,V B ) 



(1) 



where the first term is the crystal Hamiltonian of Si, the 
second and third terms the potential energy due to the 
two donor nulcei, the 4th term the nuclear repulsion of 
the two positively charged impurity cores, and the last 
term the potential energy due to the surface and barrier 
gates obtained from a commercial Poisson Solver (25|. 
The single impurity potential is expressed as, 



Voi(r) 
V Dl (r) 



-Uo, 



Ri 



r + Ri (2) 
= Ri (3) 



where Ri is the location of the i-th impurity. Interfaces 
are treated by closed boundary conditions with a model 
of surface passivation of dangling bonds (3^. The full 
atomistic Hamiltonian of about 3 million Si atoms was 
solved by a parallel Lanczos / Block Lanczos algorithm to 
obtain the ei genv alues and wavefunctions in the desired 
energy range [29[ . Each data point in this work required 
about 6 hours on 48 processors 31 1. 

The tight-binding method used here is under the hood 
of the the Nanoelectronic Modeling Tool (NEMO-3D) 
29l . 32| . and had been successfully used to verify Stark 
shift of the hyperfine coupling between a donor and its 
nucleus [l2T | in good agreement with experiments Q and 
with momentum space methods [Toj . The same method 
was also applied to investigate the orbital Stark shift of a 
donor-interface well system, and was verified with single 
donor transport experiments in FinFETs Q . 



IV. RESULTS AND DISCUSSIONS 
A. The molecular spectrum of P2+ at zero gate bias 

The lowest 6 states of a single Group V donor in Si 
are of Is type due to the six-fold degenerate conduc- 
tion band minima of Si. The sharp confining potential in 
the vicinity of the donor nucleus causes coupling between 
the valleys, a phenomena termed as valley-orbit (VO) or 
chemical splitting. The net effect of the VO interaction 
is a splitting of the 6 Is states into a singlet (Ai), triplet 
(T2) and a doublet (Ei) orbital manifolds. The strength 
of this VO interaction varies from one donor species to 
another, and is caused by a number of microscopic prop- 
erties such as variation of the dielectric constant of Si 
from its bulk value, local strain originating from bonds 
between the donor and the Si atoms, and so on. For an 



isolated P donor in Si, the (Ai), (T2) and (Ei) states are 
bound at -45.6, -33.9, and -32.6 meV respectively. When 
two hydrogenic donors are placed near each other, each 
pair of corresponding Is states gives rise to a bonding 
(symmetric) and an anti-bonding (anti-symmetric) state. 
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FIG. 2: Spectrum of the donor molecule showing the lowest 
12 molecular orbital energies as a function of donor separa- 
tion along a) [100], and b) [110]. c) and d) show the tunnel 
couplings between states 1 and 2 (A12), and 2 and 3 (A23) 
respectively. 

Fig. 2 shows the variation of the binding energies of 
the 12 lowest molecular states of the P2+ system as a 
function of impurity separation. The binding energies 
show smooth exponential decay with donor separation R 
for impurities separated along [100] (Fig. 2a), but ex- 
hibit some oscillatory behavior for impurities separated 
along [110] (Fig 2b), consistent with EMT descriptions 
[22] ], For large enough R, the bonding and anti-bonding 
pairs become almost degenerate, and the binding ener- 
gies reduce to the single P binding energies with dou- 
ble degeneracies. At small R, the bonding state arising 
from one of the T2 states of each donor approaches the 
anti-bonding state arising from Ai states, causing A23 
to be comparable in magnitude to A12. This situation 
is undesirable in quantum computer architectures since 
the qubit Hilbert space, which usually consists of the 
two lowest eigenstates, needs to be well-isolated in en- 
ergy from the rest of the manifold. A well-isolated qubit 
Hilbert space makes the qubit less prone to decoherence 
and errors during operation. In this paper, we will use 
A12 = A23 as a threshold reference level for device op- 
eration. If A12 > A23, the device is in an undesirable 
operation regime. From Fig 2a and 2b, it is observed 
that the donor separations need to be at least 7 nm for 
such coherent applications. 
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B. Surface gate control 
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FIG. 3: Surface gate response of a) A 12 and b) A23 for three 
different donor separations. The surface gate is located 10.8 
nm (Ls) from the left donor, while the donors are buried 20.1 
nm (D v ) below the oxide. 

The function of a surface gate is to provide a potential 
difference between the two impurities, and hence to con- 
trol the energy spacing of the eigen states of one impurity 
relative to the other. In the device geometry shown in 
Fig 1, a positive bias to the surface gate lowers the po- 
tential of the left impurity relative to the right impurity. 
The Ai state of the left donor moves farther below that 
of the right donor, and this increases A 12. The T2 state 
of the left donor moves closer to the Ai state of the right, 
which decreases A23 . Fig 3a and 3b show the variation of 
A12 and A23 respectively with an S-gate voltage for three 
different donor separations along [100]. For larger impu- 
rity separations, the slopes of the A i2 and A23 curves are 
steeper. If the impurity separations are larger, then the 
left impurity experiences a stronger S-gatc potential than 
the right impurity. If the impurities are closer, then the 
same surface gate voltage provides less potential drop be- 
tween the two impurities, and the surface gate response 
is weaker. 

To elucidate the different voltage regimes shown in Fig 
3, it is helpful to look at the first three eigen states of 
the whole device (top 3 rows of Fig 4) along with the 
net electrostatic potential profile the device is operating 
under (bottom row of Fig 4). Each column in Fig 4 
illustrates a snapshot of the donor molecule for a specific 
surface gate voltage with two donors separated by 15.2 
nm along [100] (points marked by square data points in 
Fig 3). 

At Vs = (col. 3), the two donor wells are aligned 
in energy, giving rise to bonding and anti-bonding pairs. 
The lowest two states (Si & S2) are formed from the 
symmetric and anti-symmetric superpositions of the Ai 
state of each donor. S3 is a bonding state arising from 
the T 2 states. 

As Vs is ramped up to 0.2 V (col. 4), the device is 
in the linear operation regime shown in Fig 3. The left 
impurity is lower in energy than the right, giving rise to 
a left localized ground state (Ai of left P) and a right 
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FIG. 4: The three lowest states of the donor molecule at 
various surface gate voltages. The bottom row shows the 
net electrostatic potential in the device, while the top 3 rows 
shows the first 3 wave function probability amplitudes. Each 
column represents the device at a specific voltage. 



localized first excited state (Ai of right P), while state 3 
arises from the T2 state of the left impurity. A plot of 
the total potential shows that the surface gate expands 
the potential contours of the left impurity relative to the 
right. 

If the S-gate voltage is reversed in polarity to -0.2 V 
(column 2), the effects described above are reversed be- 
tween the two impurities, with a right localized state ap- 
pearing as the ground state. 

At Vs = 0.5 V (col. 5), a surface well is formed near 
the gate, and quantum states begin to appear in this 
well. In a finite sized nanostructure like this, whether 
the lowest states are at the surface or at the donor is 
determined by which well is lower in energy. At Vs = 0.5 
V, the surface well is almost as deep as the donor well. 
While the ground state of the system is still seen to occur 
in the left donor, the higher states occur in the surface 
well. A small increase in bias at this point moves the 
interface well deeper than the donor well, and produces 
ionization of the donor electron. For a donor close to 
the surface, the donor well can strongly couple to the 
interface well, and give rise to the prospect of adiabatic 

ionization 0, m m, M ■ 

Since the gate confined 
surface states are closely spaced in energy compared to 
the donor states, both A12 and A23 show a sharp decrease 
at the onset of this ionization process, as shown in Fig 3. 

At Vs = —0.5 V (col. 1), the interface well is raised 
in energy, and does not play a role. The gate voltage 
is high enough to push the T 2 state of the right donor 
below the Ai of the left donor. When this happens, any 
subsequent change of the gate bias causes a small change 
in A12, and it flattens out, as shown in Fig 3a for the R = 
15.2 and 20.1 nm curves. This is expected as A12 is now 
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equal to the splitting betweeb the Ai and T2 of the the 
same donor, and has a much weaker field response. The 
A23 curve captures the energy gap between the higher 
states of the same donor, and tends to flatten out also. 
The tunnel couplings therefore saturate in this voltage 
regime. 

Since A12 is to be less than A23 for the donor device, 
the linear regimes about Vs = are important for co- 
herent operations. The operational basis can be either 
the left and right localized states at modest bias (linear 
regime) , or the bonding and anti-bonding states near zero 
bias [l|| . The ionization regime at high positive bias can 
be important if the donor electron is to be shuttled from 
one donor to another by transport along the interface 1 1511 
or if a measurement at the surface is to be performed [33| . 

C. Barrier gate control 
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FIG. 5: Barrier gate response of a) A12 and b) A23 for three 
different donor separations with the donors 20.1 nm (Dy) 
below the oxide. 

In contrast to a surface gate, a barrier gate subjects 
both donors to a similar potential. A positive barrier gate 
lowers the potential in the region between the two impu- 
rities, and increase hybridization of the left and right 
donor states. A negative barrier gate raises the poten- 
tial of this region, and provides more localization of the 
single donor states to each impurity. A12 and A23 are af- 
fected less by B-Gates than S-Gates, as exhibited by the 
smaller slopes of the Ai 2 -Vb curves of Fig 5 compared to 
those of Fig 3. The B-Gate also generates an interface 
well which eventually ionizes the donors at high enough 
gate bias, causing a sharp fall in A i2 and A23. 

Fig 6 shows how the tunnel coupling is affected when 
both S-gate and B-gate biases are present simultaneously. 
In each curve, the B-gate is held fixed at some voltage, 
and the S-gate is varied over a range of -0.5 V to 0.5 
V. It was shown in Fig 3 that at zero B-Gate bias, the 
Ai2-Vs curve shows a minimum at Vs = 0. In the pres- 
ence of a positive (negative) B-gate bias, this minimum 
shifts to a negative (positive) Vs value. Since a posi- 



tive B-gate lowers the potential in the barrier regions, 
and binds the electron more tightly between the donors, 
the symmetric-anti-symmetric gap increases. Although a 
negative S-gate bias causes detuning of the two impuri- 
ties, it also lowers the potential in the barrier region, and 
thus compensates for the positive barrier gate. Hence, 
to obtain the same zero-field symmetric-anti-symmetric 
gap, a negative S-gate bias is needed, which explains the 
shift of the Ai2-Vs curve towards a negative Vs when Vb 
is held positive. 

Since the gates are separated by 15 nm in this case, 
there is significant cross-talk between them. A positive S- 
gate therefore makes the surface well near B-gate deeper, 
and vice versa. As a result the ionization regime shifts 
to the left with increasing B-gate bias. 

Simultaneous S and B gate biases can be used to ini- 
tialize the donor molecule system. At first, each donor 
has one electron bound to it. To form a P2+ molecule, 
one electron has to removed. This can be done by hold- 
ing Vb negative so as to raise the barrier between the left 
and the right donors, while applying a simultaneous pos- 
itive bias to the S-gate to ionize the left donor electron 
to the surface well. 




Vs (V) 

FIG. 6: Effect on A12 of applying biases to the surface and 
barrier gates simultaneously. 



D. Sensitivity to Donor Placement 

In this section, we investigate the sensitivity of the tun- 
nel coupling to relative donor placements along different 
directions. In Fig 7, the radial separation between the 
donors is held fixed at 10.8 nm, while the angular sepa- 
ration is varied from to 45° at a fixed impurity depth 
of 20.1 nm. A12 and A23 are plotted in Fig 7a and 7b 
respectively for three different gate configurations. 

Although the tunnel coupling shows some oscillatory 
behavior with angular separation, both a barrier and a 
surface gate voltage are each able to change the tunnel 
coupling significantly. This shows that there is enough 
gate controllability irrespective of how the donors are 
placed relative to each other. As a note of comparison, a 
similar study was done on the two-donor exchange cou- 
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FIG. 7: a) A12 and b) A23 as a function of angular separa- 
tion between the donors. The radial donor separation and the 
donor depths are held fixed at 10.86 nm and 20.1 nm respec- 
tively, while the angle is varied from to 45° measured from 
the [100] direction. 



of Fig 8b has a steeper slope for [111] donor separations 
compared to the [100] and [110] separations. 

Similarly, for [111] donor separations, an S-gate may be 
almost equidistant from the two donors, and can act as a 
B-gate, thus reducing the voltage response of A12 signif- 
icantly. This is evident by the smaller slope of the [111] 
curve of Fig 8a compared to the other two curves. For the 
device structure considered here with donor separations 
along [111], the radial distances between the S-gate cen- 
ter and the left and right donors were 26.7 nm and 25.1 
nm respectively. This shows that donor depth below the 
gate is an important feature in experimental design as 
surface gates can act as barrier gates and vice-versa for 
certain donor orientations. 



pling J(V) in Ref |26j. In contrast to the result obtained 
here, the J(V) curve exhibited some gate controllabil- 
ity issues. For certain angular separations between the 
donors, a gate voltage was not able to alter the magnitude 
of J(V) significantly. The absence of such controllability 
issues in this case is quite encouraging. 
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FIG. 8: a) Surface and b) barrier gate response of A12 for 
donors separated along three different directions. The device 
parameters: _R=15.2 nm, Ls=10.86 nm, Dv~20.1 nm for 
R[100] and R[lll]. 

Figures 3 and 5 demonstrated smooth gate control 
over donors separated along [100]. Fig 8 demonstrates 
that such smooth control also exists for donors separated 
along other directions. In addition, Fig 8 also demon- 
strates some geometry effects that need to be considered 
in practical devices. Fig 8a shows the effect of an S-gate 
on A12 for donors separated along three different direc- 
tions, while Fig 8b shows the same for a B-Gate. 

The striking feature of the curves is the asymmetry 
in voltage control between the [111] separation and the 
other two directions. If both impurities are not on the 
same xz plane (Fig 1), which happens for [111] separa- 
tion, then a B-gate may subject one donor to a higher 
potential than the other. This is more pronounced if the 
donor depths below the gate are unequal by more than 
tens of lattice constants. In such a case, a B-gate can act 
as an S-gate, and the voltage response curve Ai2-Vb can 
show a significantly larger response. The Ai2-Vb curve 



E. Design issues: Gate placement and donor depth 
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FIG. 9: Effect of surface gate placement on A12. b) S-gate 
response of A12 for different donor depths. The impurities 
are separated by 15.2 nm along [100]. 

Fig 9a shows the effect of varying the distance between 
the S-gate and the left impurity for two donors separated 
by 15.2 nm along [100]. The slope of the A12-V5 curve is 
sensitive to L$ because it determines the potential con- 
tours the two impurities reside on. are on similar poten- 
tial contours. It is seen that the slopes of the curves are 
the steepest when Ls is between 10.8 and 20.1 nm, sug- 
gesting that there may be an optimal gate distance from 
donor for design purposes. A small L$ of or 5.5 nm 
is likely to subject both impurities to a strong potential. 
As a result, the detuning is less responsive. 

In Fig 9b, we show how the S-gate response of the 
tunnel coupling is modified for two different donor depths 
below the oxide. For both donors at a shallower depth of 
10.8 nm, the Ai 2 -"Ks curve is steeper because proximity 
of the donors to the gate lends more control. The curve 
is also seen to flatten off before reaching the ionization 
regime. This is because the stronger potential brings the 
T2 state of left impurity below Ai state of the right, and 
causes the tunnel coupling to become saturated before 
ionization is reached. Any further increase in Vs after 
0.2 V, simply captures the energy gap between the Ai 
and T2 states of the left impurity. In contrast, the curve 
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for 20.1 nm depth reaches ionization directly after the 
linear regime. 



V. CONCLUSION 

Gate control of a P2+ molecular system was investi- 
gated in detail from a tight-binding approach. It was 
found that excited states can place a limit on the range 
of operating voltages and donor separations. This arises 
from the necessity that the qubit Hilbert space needs to 
be sufficiently isolated from the other states of the system 
for robust coherent operation. Realistic TCAD [25j gate 
potentials were used to demonstrate that smooth con- 
trollability exists over the molecular states of the donors. 
The detuning of two donors by a surface gate in various 
voltage regimes was analyzed. It was also shown that the 
tunnel coupling is more sensitive to S-gatc control than 
B-gatc, with the exception of certain donor orientations 
for which the S-gate can act like a B-gate and vice versa. 
Simultaneous operation with both S and B gates was 
shown to shift the minimum energy gap to non-zero gate 
configurations. Practical control of such a system is also 
limited by interface ionization, and voltage regimes were 
established where this could take place. Sensitivity of 
the tunnel coupling to donor placement was also investi- 
gated. Although the tunnel coupling exhibits oscillatory 
behavior with donor separations in certain directions, it 
was found that a gate voltage was always able to magnify 
it significantly. This confirms that the tunnel coupling 
is free from the voltage controllability issue observed in 

The purpose of this work was to highlight the control- 
lability and design issues that an experimental imple- 
mentation and control of a P2+ molecular system would 
involve. The focus here has been not only on the physical 
trends, but also on quantitative characterization of some 
quantities that are of interest to experimentalists. This 
P2+ molecule forms the building blocks of structures that 
can utilize CTAP to transport information across qubits. 
Realization of CTAP in an actual donor array is likely to 
involve very precise gate control over quantum states of 
coupled donors. This work lays the foundation for fur- 
ther work on CTAP that involves precision gate control 
to realize an adiabatic path for population transfer across 
donor arrays. 
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